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1. Introduction 

During the past decade, there has been considerable 
interest in quantum transport phenomena in devices 
where the spin degree of freedom of the electron plays 
a crucial role e.g. the spin-transistor. 1 ) Of particular in- 
terest are mesoscopic systems where the conductance is 
sensitive not only to the sample's specific configuration 
but also to the spin population in the leads attached to 
the sample. 

One of the special characteristics of mesoscopic sys- 
tems is the fluctuating nature of transport coefficients 
such as the conductance. Such fluctuations do not de- 
pend on the details of the sample and the strength of dis- 
order but are usually determined by the symmetry and 
the dimensionality of the system. The standard classifi- 
cation 2 ) is detailed in Table. I 



Table I. Universality class and corresponding symmetry of the sys- 
tem and the properties of scattering matrix. 



Universality class 


Symmetry 


Properties of S matrix 




TRS 


SRS 




Orthogonal 


yes 


yes 


unitary symmetric 


Symplectic 


yes 


no 


unitary self-dual 


Unitary 


no 


irrelevant 


unitary 



The two relevant symmetries are time reversal sym- 
metry (TRS) and spin rotation symmetry (SRS). TRS is 
broken by an applied magnetic field or by magnetic scat- 
tering due to magnetic impurities or magnetic domain 
walls. If TRS is broken, systems are classified as unitary, 
regardless of whether or not SRS is broken. The spin- 
orbit interaction breaks SRS but preserves TRS, and in 
this case systems are classified as symplectic. 

The relation between the Hamiltonian of the sample 
and its corresponding universality class is well under- 
stood. However, the properties of mesoscopic systems 
that are measured in transport experiments depend not 
only on the sample but also on the leads through which 
currents flow in and out of the sample and through which 
voltages are measured. 

In this contribution we address the question of whether 
or not the universality class can change as a result of an 
asymmetric spin population in the leads i.e. when the 
number of up and down spin channels in the leads be- 
comes asymmetric as a result of Zeeman splitting. We 



show that even if the sample is unchanged, the univer- 
sality class can change under certain conditions. 

2. Model 

We consider a two dimensional (2D) sample connected 
to two electron reservoirs by ideal leads. The 2D sys- 
tem is constructed in the xy-plane with the x direction 
taken as the direction of current flow. The sample re- 
gion contains both a random potential and a spin-orbit 
interaction described by the Ando Hamiltonian 3 ) 



H 



E 

..j> ,a,a f 



>c\„c 



■ja' 



where 



and 



Vt 



cos 9 
— sin 9 



sin t 
cost 



V iti 



+y 



cos a i sin c 
% sin 9 cos 6 



(1) 



(2) 



(3) 



Here ct (cj CT ) denotes the creation (annihilation) oper- 
ator of an electron on site i with spin a. Wi denotes 
the random potential on the site i, distributed uniformly 
on [— W/2, W/2]. The hopping is restricted to nearest 
neighbours. The unit vector of the a;(y)-direction is de- 
noted by x(y). The parameter 9 denotes the strength of 
the spin-orbit interaction. The strength of the hopping 
is taken to be the energy unit. 

We impose the fixed boundary condition in the trans- 
verse direction y. The transverse energy of the channel 
(i, a) in the lead ef is given by 
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where n denotes the number of sites in y-direction, and 
Z denotes the strength of Zeeman splitting and a(— ±1) 
is the spin index. The channel (i, a) is a propagating 
mode if \Ep — ef\ < 2.0, Ep being the Fermi energy. For 
example, when we set Ep = —1.1 and Z = 1.0, the 
number of up (down) spin channels becomes 27 (15) for 
the sample of width 30 sites. 

In the following simulation, we set to be W — 2.0, 
9 = 0,7t/4 and Ep = —1.1. The system size is set to be 
30 x 30 in units of the lattice spacing. 
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Results 



Tabic II. Effects of Zeeman splitting in a lead 



We extend the transfer matrix method 4 - 1 to the case 
where there is a spin-orbit interaction in the sample and 
the population of up and down spins in the leads is asym- 
metric. We calculate the transmission matrix t in the two 
terminal configuration. The population of up and down 
spins in one lead is always set to be symmetric (Z = 0), 
and that in the other lead is varied (Z = 0, 1, 2). 

To detect any change of universality class, we investi- 
gate the spacing distribution P(s) where s is the inter- 
val between neighboring transmission eigenvalue r (the 
eigenvalues of tt'). An ensemble of about 10 6 samples is 
simulated for the estimation of the level statistics. 

Figure 1 shows the spacing distribution P(s) of the 
transmission eigenvalue t for the sample with a spin- 
orbit interaction. In the absence of Zeeman splitting in 
the leads (o) P(s) is close to the Wigner surmise for 
the Gaussian symplectic ensemble (GSE). On the other 
hand, with Zeeman splitting in the leads (•) P(s) is close 
to the surmise for the Gaussian unitary ensemble(GUE). 
The reason for this crossover is that the asymmetry of 
the population of up and down spins in the lead destroys 
the self-dual property of the scattering matrix. 

Figure 2 shows the spacing distribution of the trans- 
mission eigenvalue r for the sample with 6 = i.e. no 
spin-orbit interaction. In contrast to Fig. 1 the level 
statistics do not depend on the states of the leads and 
are well fitted by the formula for the Gaussian orthog- 
onal ensemble (GOE) . We summarize these results in 
Table. II. 

In summary, we have shown that the level statistics 
of the transmission coefficient is sensitive to the asym- 
metry of the spin population in the lead. We can thus 
control the conductance fluctuations without changing 
the sample parameters. For example, the variance of 
conductance is suppressed by 1/2 when we introduce the 
asymmetry of the spin population in a lead. (Note that 
the GUE in the presence of the asymmetry is spin non- 
degenerate, and the variance of conductance is 1/4 the 
spin degenerate value.) This behavior is essentially dif- 
ferent from that of the asymmetry of the number of chan- 
nels between the leads where the universality class is not 
changed. 



SO in a sample 


Zeeman in a lead 


Universality class 


yes 


no 


Symplectic 


yes 


yes 


Unitary 


no 


no 


Orthogonal 


no 


yes 


Orthogonal 



without Zeeman in a lead o 
with Zeeman in a lead • 
GSE — 
GUE - 




Fig. 1. Spacing distribution of r for the sample with the spin-orbit 
interaction. We set Z = 2. The distribution without Zeeman 
splitting in a lead (o) fits to the form of GSE while that with 
Zeeman (•) fits to GUE. 



without Zeeman in a lead o 
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Fig. 2. Spacing distribution of r for the sample without the spin- 
orbit interaction. We set Z = 1. The distribution fits GOE 
regardless of the existence of Zeeman splitting in a lead. 
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